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Abstract. We examine the predictions of the core accretion - gas capture model concerning the efficiency of planet formation 
around stars with various masses. First, we follow the evolution of gas and solids from the moment when all solids are in 
the form of small grains to the stage when most of them are in the form of planetesimals. We show that the surface density 
of the planetesimal swarm tends to be higher around less massive stars. Then, we derive the minimum surface density of 
the planetesimal swarm required for the formation of a giant planet both in a numerical and in an approximate analytical 
approach. We combine these results by calculating a set of representative disk models characterized by different masses, sizes, 
and metallicities, and by estimating their capability of forming giant planets. Our results show that the set of protoplanetary disks 
capable of giant planet formation is larger for less massive stars. Provided that the distribution of initial disk parameters does 
not depend too strongly on the mass of the central star, we predict that the percentage of stars with giant planets should increase 
with decreasing stellar mass. Furthermore, we identify the radial redistribution of solids during the formation of planetesimal 
swarms as the key element in explaining these effects. 



1. Introduction 

Radial velocity surveys led to the discovery of over 150 extra- 
solar planets around main sequence stars. Published descrip- 
tions of most of them ca n be found in the r eferences given 
bv lMarcv et al.l J2005t) and lMavor et alJ J200 4). Those surveys 
have been the most successful in the case of G dwarf stars, be- 
cause such stars have well-defined spectroscopic features and 
show only a little photospheric activity. Consequently, most of 
the known extrasolar planets orbit stars similar to our Sun. Due 
to the constant progress in the detection techniques, the ob- 
servational programs recently started to also include stars with 
lower masses on a larger scale, namely M dwarfs. Moreover, 
some of these s urveys are now particularly dedicated to lower- 
mass stars (e.g. lEndl etaD 120031: lionfils et al.ll2004 . So far, 
these efforts have led to the disc overy of three planets ar ound 
two M d warf stars: Gliese 8 76b,c Jlvlarcv et al.ll998Ll200lh and 
GJ436b feufler et all2004 . 

From the theoretical point of view, the problem of giant 
plan et formation arou nd M dwarfs was studied recently by 
lLaughlin et alJ d2004h . They addressed it within the core ac- 
cretion - gas capture model (CAGCM) that provides the most 
widely accepted scenario explaining the formation of giant 
planets in both the Solar System and extrasolar planetary sys- 
tems. This model predicts that first a solid planetary core is 
formed by collisional accumulation of planetesimals. When the 
core reaches a mass of a few Earth masses, it starts to accrete 
gas, and an extended hydrostatic envelope is built around it. 
As the accretion rate of gas is greater than the accretion rate 
of solids at this time, the envelope eventually becomes more 



massive than the core. When this happens, a runaway accre- 
tion of gas ensues, which is terminated either by tidal interac- 
tions of the planet with the protoplanetary disk or by the dis- 
sipation of the disk. CAGCM has found supporting evidence 
in the discovery that s tars with planets have higher metal- 
licities than field stars ( Sa ntos et al.ll2000t iFischer & Va lenti 
2003j). This is because the formation time of giant planets 
decrease s with increasing s urface density of the planetesimal 
swarm (Polla ck et alJ ll996). which in turn increases with the 
primordial metallicity of the protoplanetary disk. Thus, giant 
planets are expected to form more easily in disks with higher 
metallicity. 

lLaughlin et alJ J2004 conclude that M dwarfs have a lim- 
ited ability to form Jupiter-mass planets. This is a direct conse- 
quence of their assumption that the surface density of the plan- 
etesimal swarm out of which planetary cores are formed scales 
linearly with the mass of the central star. However, the solid 
component of the protoplanetary dis k evolves in a different way 
than the gaseous component JWeidenschilling & Cuzzil [T993). 
Due to the gas drag, a significant redistribution of solids takes 
place, and in the inner disk their surface density can be substan- 
tially enhanced compared to the i nitial one JWeiden schilling 
2003; Stepinski & Valageas l!997h . In general, the efficiency 
of the processes responsible for the redistribution of dust de- 
pends on the mass of the central star. An obvious conclusion 
is that analysis of the formation of giant planets around stars 
with various masses should include the global evolution of 
solids in protoplanetary disks. A simple m odel o f fhe evolution 
of so lids was proposed by iKornet et alJ J2002I) : iKornet et alJ 
(2005). Applying it to solar-like central stars, these authors re- 
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produced the observed correlation between stellar metallicity 
and the probability of a pl anet occurring (a similar result was 
independently obtained bv llda & Linl200 4b). 

The rapid progress in observational techniques opens up 
the possibility of testing the correctness and predictive power 
of the model proposed by Kornet et al To that end, we 

extend their analysis and calculate probabilities of planet oc- 
currence around stars with different masses (both smaller and 
larger than 1M Q ), which may be compared to future observa- 
tional data. In Sect.|2]we explain our approach to the evolution 
of protoplanetary disks and planet formation. The results of our 
calculations are presented in Sect.[3]and discussed in Sect.@] 

2. Methods of calculation 

2.1. The disk 

We model the protoplanetary disk as a two-component fluid, 
consisting of gas and solids. The ga seous com ponent is de- 
scribed by the analytical model of Ste ninskil lll998l) . which 
gives the surface density of gas, S g , as a function of distance 
a from the star and time t, in terms of a selfsimilar solution 
to the viscous diffusion equation. The viscosity coefficient is 
given by the standard a prescription: 



(1) 



where H is a density scale-height of gaseous disk and Cs de- 
notes the speed of sound in the gas. All other quantities charac- 
terizing the gas are obtained in a thin disk and vertical thermal 
balance appro ximation by solving the set of equations (see, for 
example Fran k et all 19921) : 





= 2H Pg 


(2) 


H 


V2C S 
Q K 


(3) 


cl 


k B T 
lum H 


(4) 


16cr B r 4 

3kL„ 


= -Z„v£ll . 

4 = K 


(5) 



Here T is the temperature in the midplane of the disk, and Q.^ 
the Keplerian angular velocity. For the Rosseland mean opacity 
k the analytical piecewis e-continuous power law formulas from 
iRuden & Pollackl Jl99ll) were adopted. In the parts of the disk 
wher e the optical thickness fa lls below a critical value of T CI j t = 
1.78 IRuden & P ollack 1991), the last equation is replaced by: 
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The main assumptions underlying our approach to the evo- 
lution of solids are the folllowing. (1) At each radial distance 
from the central star the particles have the same size (which 
in general varies over time). (2) There is only one compo- 
nent of dust, in this case corresponding to high-temperature 
silicates with the evaporation temperature r eva p = 1350 K and 
bulk density 3.3 g crrT 3 . This choice is justified by the fact 



that in most cases the surface density of a water-ice planetesi- 
mal swarm would be too low in the range of distances consid- 
ered here, e.g. r < 5 AU, to e nable formation of a giant planet 
in a reasonable time (see also Kornet et al. 2004). We include 
other species, e.g. water ice only, in modeling the gaseous disk 
as the sources of opacity. In other words we assume that the 
evolution of silicate grains is not strongly influenced by the 
evolution of i.e. ice grains. (3) All collisions between particles 
lead to coagulation. (4) When the temperature exceeds r evap , 
local solids immediately sublimate and the vapour evolves at 
the same radial velocity as the gas component. (5) Initially, 
in regions where the disk temperature is lower than r evap , all 
solids are in the form of grains with radii 10~ 3 cm (the results 
do not depend on the choice of that particular value, as long 
as the solids are initially small enough to be coupled well to 
the gas). (6) The radial velocities of solid particles are entirely 
determined by the effects of gas drag. (7) The relative veloci- 
ties of solid particles when they collide are computed accord- 
ing to the turbulent model described by Stepin ski & Valag eas 
(8) The evolution of the solids does not affect the evo- 
lution of the gas. At each radius, the vertical extent of the solid 
particle distribution is calculated and is evolved in time, so the 
effect of the sedimentation of solids toward the midplane of 
the disk is taken into ac count. All assumption s and approxima- 
tions a re di scussed in IStepinski & Valaeeasl dl997h . IStepinskil 
d!998l) . and lKornet et alJfcOOlb . 

The evolution of solids is governed by two equations. The 
first of them is the continuity equation for the surface density 
of solid material, E s . The second one, describing the evolution 
of grain sizes, can be interpreted as the continuity equation for 
size-weighted surface density of solids, Z s i ze = ,s(a)E s , where 
s(a) is the radius of solid particles at a distance a from the star. 
The equations are solved numerically on a moving grid whose 
outer edge follows the outer e dge of the dust dis k. The details 
of the method can be found in lKornet et al.l d200ll) . 

2.2. The planets 

We model the formation of a giant planet in situ, so the orbital 
parameters of the planet do not vary in time. Our procedure for 
the evolution of the protoplanetary cores is based on the follow- 
ing assumptions: (1) core accretion starts when solids at a given 
distance a from the star reach radii of 2 km; (2) at each time, the 
planetesimals are mixed well through the feeding zone of the 
planet, so their surface density 2 S is always uniform in space, 
but usually decreasing with time as planetesimals accrete onto 
the planet; (3) the planetesimals do not migrate into the feed- 
ing zone from outside or inside and vice versa, but they can 
be overtaken by the boundary of the feeding zone as it expands 
due to the growing mass of the planet. Under these assumptions 
the growth of protopl anetary core mass M c can be d escribed by 
the formula given by Papaloizou & Terauem ( IS 
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is the radius of the Hill sphere of the planet. The quantity h 
is a constant factor that reflects different definitions of the Hill 
radius in the literature. Herein we assume h - 3. The value of 
C\ given bv lPapaloizou & TerauerrJ (^999) is 81^/32; we use 
a factor of 5 ( the difference comes from a different definition 
of Ru). The quantity C cap describes the increase in the effective 
capture radius of the planet with respect to its real radius R p due 
to interaction o f planetesimals with the envelope of the planet 
( Podolak et al. 1988). We approximate it with a fit to the re sults 
of Bodenheimer et al. ( 2000) provided by Hubickvi ( 2001 ). For 
core masses less than 5M , no increase in the effective capture 
radius is assumed, i.e. C cap (M c < 5M e ) = 1. For higher core 
masses, it increases linearly with the mass of the core, reaching 
its maximum value of C cap = 5 for M c = 15M e . The surface 
density of planetesimals E s also changes in time, as they are 
accreted onto the core, and the feeding zone expands. Under 
our assumptions it can be calculated as 
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where £ Si i n j t is the initial surface density of planetesimals, and 
Aa is the width of the fe eding zone. Herein we assume that 
Aa = 8# H dLissaueJl993l) . 

To calculate the rate of gas accretion onto the protoplanet, 
one should solve the equations of mass conservation, hydro- 
static equilibrium, energy generation from accretion of plan- 
etesimals and quasi-static contraction, and radia tive or convec- 
tive en ergy transp ort, as given e.g. in iBodenheimer & Pollack! 
( 1986). Following Id a & Linl J2004ah . we use a simplified ap- 
proach based on fits to the numerical results. We assume that 
the accretion of gas starts when the core reaches critical mass 
of 
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Jlkoma et all20 00). where K em is the opacity in the envelope of 
the planet. Its actual magnitude is currently poorly constrained. 
We assume that K env = 1 cm 2 g -1 . For the value of the powerlaw 
index in the dependence on M c (Eq.I10> we adopt 0.25. When 
the mass of the protoplanet is higher than M c cr j t it contracts 
on the Kelvin-Helmholtz ti me scale t kh- Brvd en et alJ (12000) 
show by fitting the result of Pollac k et al.l 11199 6) that 
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where b — 10 and c = -4. Consequently, for M p > M ccl -i t we 
adopt the following equation for the gas accretion rate onto the 
planet 



dM env M p /M p . . 
— — = — = A — M e yr 1 = am: 
At T KH I M ffl 1 
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where A = 10~ 10 . The free parameters in Eqs. (I7l- I12> are the 
mass of the central star M*, the distance a of protoplanet from 
the star, and the initial surface density of planetesimals S s> i n i t . 



3. Results 

In the CAGCM the process of planet formation is naturally split 
into two main phases. In the first one, dominated by collisional 
accumulation of dust grains, a planetesimal swarm is formed 
in the protoplanetary disk. In the second phase, dominated by 
gravitational interactions, planetary cores are assembled and 
subsequently accrete planetesimals and gas from the disk. In 
the following subsections we investigate how each phase is in- 
fluenced by the mass of the central star. 



3.1. From dust grains to planetesimals 

To illustrate how the mass of the central star influences the for- 
mation of the planetesimal swarm, we follow the evolution of 
a protoplanetary disk with an initial mass Mq = 0.1 M and 
initial outer radius R Q = 1000 AU for three values of M* (0.5,1 
and 4 M ). Figure [l] shows the initial distributions of E g and 
T in the midplane of the disk. In all three cases the gas is dis- 
tributed very similarly, with E g dropping monotonically from 
~ 4000 g ctrT 2 at 0.01 AU from the star to 0.1 g ctrT 2 at the 
outer edge of the disk. The changes of slopes in the distribution 
of Eg correspond to transitions between different powerlaws de- 
scribin g the opacity of th e disk matter in different temperature 
ranges ( lRuden&Pollacklll99lh . The distribution of tempera- 
ture is qualitatively very similar to the distribution of £ g . In the 
disk around a 1 M Q star T drops from 2500 K at 0.01 AU to a 
few Kelvins at the outer edge of the disk. The evaporation tem- 
perature assumed for the dust grains in our models (1350 K) is 
reached at 0.05 AU. Note that at a given radius T increases as 
the mass of the central star is increased. It is the result of the 
increasing vertical component stellar gravity, due to which the 
scale height of the disk is reduced. 

Initially, the dust is well mixed with the gas, with the ra- 
tio 2 s /2 g = 6 x 10~ 3 being constant everywhere in the disk. 
As the disk evolves, surface density and temperature of gas 
slowly decrease due to accretion and viscous spreading; how- 
ever, the dust component evolves in quite a different way. The 
grains grow in size due to mutual collisions and gain inward 
radial velocities due to the gas drag. If they cross the evapora- 
tion radius, they sublimate and are accreted onto the star on the 
viscous timescale as a vapour. However, if their growth time 
is shorter than the timescale of inward migration, they man- 
age to reach sizes of a few km before reaching the evaporation 
radius. Their radial motions are then stopped and the planetesi- 
mal swarm attains its final form. Figure[2]shows the distribution 
of planetesimals in our model after 10 6 yr from the beginning 
of its evolution. In all three cases, the outer radius of the plan- 
etesimal swarm is much smaller than the initial outer radius of 
the disk. The difference is larger for models with smaller M*. 
This is because small solid bodies evolving in disks around less 
massive stars gain higher inward velocities and tend to migrate 
to smaller radius before reaching km-sizes. The maximum ve- 
locity of the inward drift can be estimated as 
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Fig. 1. The initial surface density (left panel) and temperature (right panel) as a function of distance from the star in the proto- 
planetary disk with an initial mass Mq =0.1 M e and outer radius Rq = 1000 AU. Different curves correspond to different values 
of the mass of the central star in solar masses as described by the labels. 
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Fig. 2. The surface density of solids as a function distance from 
the star in the protoplanetary disk with the same parameters 
as on Fig. Rafter 10 6 yr from the beginning of its evolution. 
Different curves correspond to different values of the mass of 
the central star in solar masses as described by the labels. 

JWeidenschilhngll 1 977t iKornet et al.ll2002l) . A change M* in- 
fluences it in two opposite manners. On the one hand, V™ ax 
increases with the mass of the star due to the increase in the 
disk temperature . On the other hand, it decreases because of 
the increase in the Keplerian velocity. Of these two competing 
factors, the second one dominates, and as a result the V™ ax is a 
decreasing function of M*. 

Due to the inward migration of solids and their confinement 
to much smaller radii, the final surface density of planetesimals 
is increased locally within a factor of a few in comparison with 
the initial value of S s . As this effect is larger in disks around 
less massive stars, their final planetesimal swarms tend to be 
more favourable to the formation of giant planets. 

3.2. Minimum surface density 

To quantify the influence of M+ on the effectiveness of giant 
planet formation from a planetesimal swarm, we introduce the 



concept of the minimum surface density £ s , m ; n . We define it as 
the minimum value of the initial surface density of planetesi- 
mals S s i n i t needed to form a Jupiter-mass (lMj) planet in less 
than the lifetime of the protoplanetary disk Tf. For Tf we adopt 
a value of 3 x 10 6 yr. 

First, by solving the set of equations @, 0, and H2\ 
with different values of S s i n i t , we determine £ sm j n as a func- 
tion of distance from the star. The results are shown in Fig. 
13 for the same values of M* as before. Close to the star 
(a < 10AU), £ s ,init is a decreasing function a. In this regime, the 
accreting protoplanetary core rapidly accumulates all planetes- 
imals in its feeding zone and reaches the isolation mass M; so . 
Afterwards, the accretion of planetesimals is negligible and the 
planet grows mainly due to the accretion of gas. Upon integrat- 
ing Eq. dl 2L we obtain the minimum isolation mass needed to 
form a lMj planet within Tf: 



M iso > [M]~ c + (c- 1)At { ] 



(14) 



where we assumed that the time needed for the core to reach 
Mjsc is much shorter than Tf. On the other hand, from the defi- 
nition of M; so , we see that 
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Combining these two equations and neglecting the factor Mj , 
we obtain an analytical formula for S s>m i n 



2n 1 
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(16) 



(I16> which for a < 10 AU agrees very well with the exact 
values of S s m i n on the descending branch of the curve £ Sitn j n (a) 
(see Fig. |4j. Note that in this regime the quantity E s , m ; n is an 
increasing function of the stellar mass. This behaviour is due to 
the fact that, while we increase M*, the Hill sphere shrinks, and 
as a result the isolation mass decreases (provided, of course, 
that 2 Sj i n it stays constant). 

For sufficiently large radii (a > 10 AU), E S;m i n (a) changes 
its slope and becomes an increasing function of a. In this 



K. Komet et al.: Formation of giant planets around stars with various masses 



5 



1000 



100 r 




100 



Fig. 3. The minimum surface density of a planetesimal swarm 
needed to from a lMj planet in less than 3 x 10 6 yr as a function 
of distance from the central star. Different curves are obtained 
for different masses of the central star, as labeled in units of 
solar masses in the upper right corner. 

regime, the time scale of accretion of the planetesimals onto 
the core is larger than the lifetime of the disk and becomes 
the main factor to determine £ s ,min- Consequently, the surface 
density of planetesimals in the feeding zone never drops much 
below its initial value. To describe the formation of a planet an- 
alytically under these conditions, we divide the whole process 
into two phases. During the first phase the planet exclusively 
due to the accretion of planetesimals grows . We assume that 
the surface density of planetesimal swarm S s does not change 
in time, because the core only accumulates a negligible fraction 
of solids present in the feeding zone. From Eq. we obtain 



dMp _ 
dr 
where 



A C M: 
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M 
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is a constant factor (we neglect the changes of C cap in time, 
and p c in the above expression denotes the density of the core). 
During the second phase the planet grows exclusively due to 
the accretion of gas at a rate described by Eq. dl 2i . This phase 
begins when the planet reaches the mass M p ; nt for which the 
accretion rate of gas is equal to the accretion rate of planetesi- 
mals. From Eqs. (I12> and fl!7i . we get 



M, 



p,int 
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3/(3c-2) 



(19) 



Integrating Eqs. and d!2t over time yields the lengths of 
both phases At\ and A?2: 
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Fig. 4. The minimum surface density of the planetesimal swarm 
needed to form a lMj planet in less than 3 x 10 6 yr as a func- 
tion of distance from the central star, calculated with a constant 
value of C cap = 5. The solid line represents the results of exact 
values obtained numerically, while dotted lines represent ana- 
lytical approximations given by Eqs. d!6l > and J22l >. The mass 
of the central star is 1 M . 



Combining Eqs. ( 118ft . ( 119b . 120t . and ( 12 1 i with the condition 
that Afj + Af2 < Tf , we get the equation for £ s , m j n in that regime: 
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(22) 



With Ccap = 1-2, the above formula satisfactorily reproduces 



the rising branch of curve 2 s m i n (a) (see Fig. 0}. Note that the 
"best-fit" value of C cap is consistent with the maximum and 
minimum values of this parameter used in our numerical cal- 
culations (1 and 5, respectively). At large distances from the 
star, £ s , m i n is a decreasing function of M*, because increasing 
the mass of the star causes the accretion rate of planetesimals 
to decrease (see Ea.ll7>. 

As we see, the minimum surface density of the planetesimal 
swarm required for the formation of a gas giant planet in a time 
shorter than the lifetime of protoplanetary disk is a complicated 
function of M* . In a given planetesimal swarm such a planet 
forms more easily around less massive star if its orbital radius 
is smaller than ~ 10 AU, and around a more massive star if its 
orbital radius is larger than ~ 10 AU. 

3.3. Grid of models 

The results of the last two subsections allow us to investigate 
the influence of M* on the whole process of giant planet forma- 
tion for a broad set of models of protoplanetary disks. We cal- 
culate the grid of models similar to the one described in Sect. 
13.11 but with different values of the initial disk mass Mo and 
outer radius Ro. To cover the range of masses and sizes of disks 
observed in nature, we choose Mo in the range of 0.02 to 0.2 
M Q . The range of Ro is adjusted for every metallicity so that all 
models in which formation of giant planets is possible could be 
accounte d for. For the viscosity coeffi cient a we chose a value 
of 0.001 (Panaloiz ou & Nelsonll20 03). We follow each model 
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Fig. 5. The plane of initial parameters of protoplanetary disk 
models [Mo, log Rq]. Minimum and maximum distance from 
a 1M star at which the formation of a lMj planet is possible 
within 3 x 10 6 yr is indicated by contours and a grey scale, 
respectively. White circles indicate disk models in which the 
surface density of the planetesimal swarm is everywhere lower 
than the critical value for planet formation. Black circles indi- 
cate disks in which all solids are accreted onto the star. 



until all solids are in the form of planetesimals or are accreted 
onto the star. Then, we evaluate every model with planetesi- 
mals to determine whether the surface density of planetesimals 
exceeds the minimum surface density 2 Sjm i n anywhere in the 
disk. Models with this property are labeled as planet bearing. 
For each such model we determine the minimum and maxi- 
mum distances from the star at which the surface density of the 
planetesimal swarm is larger than E,^ . The results obtained 
for different values of the stellar mass are shown in Fig. 15171 

The area of the region occupied on the [Mo, Rq] plane by the 
planet bearing disks is clearly anticorrelated with the mass of 
the central star. As we have shown in Sect. 13.11 solid grains gain 
higher inward velocities in disks around less massive stars, and 
the resulting planetesimal swarms have higher surface densi- 
ties. We see that reducing the mass of the central star increases 
the maximum Rq for which planet formation is possible in disks 
with the same initial mass Mq. 

We can also see that in a given disk the inner edge of the 
planet-bearing region moves inward as we decrease the mass 
of the star (for M+ = O.5M its radius is on the average equal 
to ~ 0.55 that for M* = 1M Q ). This effect is mainly caused 
by the larger surface density of planetesimal swarms produced 
by disks around less massive stars. Differences in X s , m i n , while 
appreciable, are much less important (when M* is reduced by 
a factor of 2, £ Sjln i n drops by 20% only; see Ea.ll6>. The outer 
edge of the planet-bearing region also moves closer to the star 
(for M+ = O.5M its radius is on the average equal to ~ 0.75 
that for M* = 1M ). This is because in most cases it coincides 
with the outer edge of the planetesimal swarm, which is more 
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Fig. 6. Same as Fig.|5]but for the central star with a mass of 
O.5M . 
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Fig. 7. Same as Fig.|5]but for the central star with a mass of 
4M Q . 



compact around less massive stars. The difference in the mini- 
mum surface density E s ,min also tends to decrease the outer ra- 
dius of the planet-bearing region, but it is again a second-order 
factor. 

Generally, our model predicts that giant planets tend to 
form at tighter orbits around less massive stars, and wider or- 
bits around more massive stars. However, at least in some cases 
their locations may be influenced by the effects of migration. 
We return to this point in Sect. |4] 
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3.4. Metallicity relation 

One of the main results of extrasolar planet searches is the 
discovery that planet-be aring stars tend to have higher metal- 
licitie s than field stars dSantos et alJl200ot iFischer & Valentil 
2003]). That correlation can be easily explained within 
CAGCM. In this scenario, the formation time of giant planets 
decreases with increasing surface density of the planetesimal 
swarm (see Eqs. (I16> and d22» . which in turn is an increas- 
ing function of the original metal content of the protoplanetary 
disk. Consequently, the giant pla nets form mor e easily in disks 
with higher metallicities. iKornet et alJ d2005l) calculated the 
rates of giant-planet occurrence in disks with different metal- 
licities around stars with mass M* = 1M . Their approach to 
the evolution of solids and to the formation of giant planets 
was the same as the one used in this paper. They were able to 
reproduce the observational correlation for disk models with 
viscosity parameter a = 10~ 3 - 10~ 2 . Herein we extend their 
calculations onto disks around stars with various masses. 

The change in the disk metallicity influences the processes 
leading to the formation of planets in two ways. First, it 
changes the structure of the gaseous disk by changing the opac- 
ity in the disk. In our models we scaled the opacity by a con- 
stant factor Z equal to the metallicity of the disk expressed in 
solar units. This approach is justified by the fact that the opac- 
ity in protoplanetary disks is mainly due to dust grains and 
molecules. Second, the primordial metallicity of the disk de- 
termines the initial ratio of dust-to-gas surface densities. In our 
models this ratio is initially independent of the distance from 
the star and is equal to 



2 s /2 e = 6 x 1(T 3 Z. 



(23) 



We compute grids of models similar to those described in 
Sect. l3.3l for eight different values of Z distributed between 0.2 
and 3. Additionally, we check the gravitational stability of the 
corresponding gaseous disk for every value of [Mo,Rq] . In 
some cases the value of the Toomre parameter 



Q 



ttGI ' 



(24) 



where S = 2 g +2 S , drops below 1 in the outer region of the disk, 
which means that they are unstable with respect to axisymmet- 
ric modes. We assume that such region fragments and giant 
planets are formed there on a very short time scale, consuming 
and/or dispersing the unstable part of the disk. In such cases we 
use modified initial values of Mq and Rq, which correspond to 
the mass and outer radius of the stable part of the original disk. 
In principle, some "special treatment" should also be applied 
to disks with 1 < g < 1.3, which develop spiral arms and for 
a while evolve due to gravitatio nal rather than viscous angu - 
lar momentum transfer (see e.g. Laughlin & R ozvczkal ll996). 
However, since the efficiency of gravitational transport quickly 
decreases with the ratio Mq/M* and it might become signifi- 
cant in our study only for the most massive disks around 0.5 
M Q stars, we decided to neglect this effect altoget her. 

Following the procedure described by Kor net et alJ d2C 
for every Z we calculate the area A p> 5 of the region occupied on 
the [Mo,log/?n] plane by disks that form planets at distances 




Fig. 8. The rate of planet occurrence as a function of the pri- 
mordial metallicity of protoplanetary disks. Different lines are 
obtained for models with different masses of the central star, as 
labeled in solar units in the upper left corner. The histogram 
shows the observational data compiled by IFischer & Valenti 
ll2003h . 



smaller than 5 AU from the central star. The last restriction 
comes from the fact that currently we k now only one ext raso- 
lar planet on a larger orbit 1 - 55 Cnc d JMarcv et al.ll2002l) . A 
measure of the rate of planet occurrence can be defined as 



P -^1 
p " C ' 



(25) 



where the normalization factor C is chosen in such a way as to 
reproduce the observed value of P p for Z = 3 and M* = 1 M Q . 

The results are presented in Fig. [S] It shows the rate of 
planet occurrence as a function of disk metallicity for three val- 
ues of Mi,. As expected, P p is an increasing function of Z. The 
minimum value of Z below which no giant planets are formed 
at orbits smaller than 5 AU decreases with the mass of the 
central from Z m ; n * 0.6 for M* = O.5M , to Z m j n m 0.2 for 
M+ = 4M . 

We see that for Z smaller than ~0.2, P p is an increasing 
function of M*. This is because most disks in which planets 
formation would be possible around less massive stars have 
outer parts that are gravitationally unstable, and the amount 
of solids present in their stable parts is too low to enable 
subsequent formation of giant planets according to CAGCM. 
However, as Z increases, smaller and smaller disks become 
planet-bearing for every Mo, and the percentage of gravita- 
tionally unstable disks in which formation of giant planets is 
not possible decreases. Consequently, the factors promoting the 
formation of a giant planet around less massive stars as de- 
scribed in previous sections become important, and P p changes 
into a decreasing function of M*. 



1 see The Extrasolar Planets Encyclopedia at 
http://www.obspm.fr/encycl/encycl.html 
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4. Conclusions 

Based on a simple approach to the evolution of solids in pro- 
toplanetary disks, we investigated the influence of the mass of 
the central star on the formation of giant planets. We showed 
that due to the more efficient redistribution of solids the plan- 
etesimal swarms around less massive stars tend to have higher 
surface densities. Next, we derived the minimum surface den- 
sity of the planetesimal swarm needed to enable formation of a 
giant planet within the lifetime of the protoplanetary disk, and 
we found that at distances from the star smaller than ~ 10 AU it 
increases with the stellar mass. Farther away from the star the 
minimum density becomes anticorrelated with the mass of the 
star. However this effect is offset by the anticorrelation men- 
tioned already between the mass of the star and the surface 
density of the planetesimals. 

These two effects determine the set of initial parameters 
characterising protoplanetary disks that are capable of giant 
planet formation within the core accretion - gas capture sce- 
nario. We showed that this set is larger for less massive stars. 
This means that the percentage of stars with massive planets 
should increase with decreasing stellar mass (at least in the 
range 0.5 M G - 4 M ). However, as discussed below, in the 
currently accessible range of orbital radii (< 5 AU), the situa- 
tion is not all that clear. 

Based on the sets obtained for different metallicities, we 
determined the occurrence rate of planets with orbits smaller 
than 5 AU as a function of the mass and metallicity of the 
star. We took into account the fact that the outer region of the 
disk is gravitationally unstable in some models. Such regions 
are located farther than 5 AU from the central star, and planets 
formed there by disk fragmentation are not included in our oc- 
currence rate. However, their presence reduces the amount of 
solid material available for the formation of planetesimals. For 
less massive stars this effect is so strong that it overcomes fac- 
tors promoting planet formation, so that for metal-poor disks 
the rate of planet occurrence decreases with the mass of the 
central star. As a result, the minimum metallicity at which gi- 
ant planets can form at orbits smaller than 5 AU decreases from 
~ 0.6 for stars with masses of O.5M to ~ 0.2 for 4M . 

In the metal-rich regime the percentage of entirely stable 
disks in which formation of giant planets is possible is larger, 
and stable regions of partly unstable disks contain enough 
solids to produce planetesimal swarms capable of giant planet 
formation. Consequently, both factors promoting planet forma- 
tion around less massive stars are in play, and a clear anticor- 
relation between the stellar mass and planet occurrence rate is 
observed. At the same time, our model does not account for the 
presence of giant planets around metal-poor stars. This may be 
due to the fact that we do not include planets that have formed 
beyond 5 AU and later migrated inward. Such an assumption 
is valid as long as the number of these planets is small com- 
pared to the number of planets that have formed within 5 AU. 
However, as we move to lower metallicities, the percentage of 
giant planets with silicate cores decreases (the silicates sim- 
ply become too scarce), while the percentage of planets form- 
ing from ice grains increases. Thus, in metal-poor systems the 
number of planets with ice cores that migrated from large orbits 



can become a large fraction of planets at orbits smaller than 5 
AU. 

Obviously, our description of the evolution of solids is very 
simplified. The basic underlying assumptions like the single- 
size distribution of solid grains or the n eglect of pla n et mi- 
grati on already have be en discussed by iKornet etal ] d2004h 
and lKorneteTail §005). The main additional assumption in- 
troduced in the present paper is the independence of the initial 
parameters of protoplanetary disks on the mass of the central 
star. While admitt edly ad hoc, i t seem s to be better than the 
one adopted by Laug hlin et alJ ( 120041) . who scaled their ini- 
tial surface density of planetesimals linearly with the mass of 
the star. They did not take into account the antecedent evolu- 
tion of solids leading to the formation of planetesimal swarms, 
and concluded that the formation of giant planets around low- 
mass stars is difficult. Recent observations suggest that masses 
of protoplanetar y disks do not st rongly depend on masses of 
the central stars (Guilloteau 2005). Nevertheless, to investigate 
the influence of our assumption, we performed additional set 
of calculations with a mass of the central star of 0.5M o and 
with the initial masses of disks scaled by factor of 0.5. The re- 
sults are shown in Fig. [9] In this case the probability of finding 
a planet does not seem to depend strongly on the mass of the 
central star, which is true for the whole range of metallicities 
we have considered. Still, our models show that the evolution 
of solids leading to the formation of planetesimal swarms is a 
vital factor facilitating the formation of giant planets, whose 
role should be particularly clear for low-mass stars. 

Our mode ls of gaseous disks do not reproduce recent ob- 
servations by (Muzerol le et alJl2005l) . which show that the ac- 
cretion rate in protoplanetary disks increases with the mass of 
the central star. However, in the mass range considered here 
this dependence is very weak, and for a given value of stellar 
mass the spread in accretion rates reaches two orders of mag- 
nitude. In our opinion these data do not invalidate our basic 
assumption that initial disk parameters do not depend on the 
mass of the star. We also assumed that heating by stellar ra- 
diation is negligible, whereas at least in some cases it can be 
a dominant source of energy in the outer regions of the disk 
(more efficient than the turbulent dissipation). As such, it may 
substantially change the structure of the disk and the radial ve- 
locities of solids. Currently we are working on models that will 
take these effects into account. 

Acknowledgements. This project was supported by the German 
Research Foundation (DFG) through the Emmy Noether grant WO 
857/2-1 and the European Community's Human Potential Programme 
through the contract HPRN-CT-2002-00308, PLANETS. KK and MR 
acknowledge support from the grant No. 1 P03D 026 26 from the 
Polish Ministry of Science. 

References 

Bodenheimer, R, Hubickyj, O., & Lissauer, J. J. 2000, Icarus, 
143, 2 

Bodenheimer, P. & Pollack, J. B. 1986, Icarus, 67, 391 
Bonfils, X., Delfosse, X., Udry, S., et al. 2004, in Astronomical 

Society of the Pacific Conference Series, 101 
Bryden, G., Lin, D. N. C, & Ida, S. 2000, ApJ, 544, 481 



K. Kornet et al.: Formation of giant planets around stars with various masses 



9 



0.3 
0.25 
0.2 - 



observations 
M=1 
M=0.5 
M=0.5, m scaled 




Fig. 9. The rate of planet occurrence as a function of the pri- 
mordial metallicity of protoplanetary disks. Disks around stars 
with masses of 1M are represented by a solid line. Disks 
around stars with masses of 0.5M o are represented by dotted 
and dashed lines. In the first case the range of initial masses of 
the disks is the same as in the case of solar type stars, while in 
the second it was scaled according to the mass of the central 
star. 



Marcy, G. W., Butler, R. P., Vogt, S. S., et al. 2005, ApJ, 619, 
570 

Mayor, M., Udry, S., Naef, D., et al. 2004, A&A, 415, 391 
Muzerolle, J., Luhman, K. L., Briceno, C, Hartmann, L., & 

Calvet, N. 2005, ApJ, 625, 906 
Papaloizou, J. C. B. & Nelson, R. P. 2003, MNRAS, 339, 983 
Papaloizou, J. C. B. & Terquem, C. 1999, ApJ, 521, 823 
Podolak, M., Pollack, J. B., & Reynolds, R. T. 1988, Icarus, 73, 

163 

Pollack, J. B., Hubickyj, O., Bodenheimer, P., et al. 1996, 

Icarus, 124, 62 
Ruden, S. P. & Pollack, J. B. 1991, ApJ, 375, 740 
Santos, N. C, Israelian, G., & Mayor, M. 2000, A&A, 363, 228 
Stepinski, T. F. 1998, Icarus, 132, 100 
Stepinski, T. F. & Valageas, P. 1997, A&A, 319, 1007 
Weidenschilling, S. J. 1977, MNRAS, 180, 57 
Weidenschilling, S. J. 2003, in Lunar and Planetary Science 

Conference #34, abstract #1707 
Weidenschilling, S. J. & Cuzzi, J. N. 1993, in Protostars and 

Planets III, ed. E. Levy & J. Lunine (Tuscon: Univ. of 

Arizona Press), 1031 



Butler, R. P., Vogt, S. S., Marcy, G. W., et al. 2004, ApJ, 617, 
580 

Endl, M., Cochran, W. D., Tull, R. G, & MacQueen, P. J. 2003, 
AJ, 126, 3099 

Fischer, D. & Valenti, J. A. 2003, in Scientific Frontiers in 
Research on Extrasolar Planets, ASP Conference Series, Vol 
294, ed. D. Deming & S. Seager (San Francisco: ASP), 1 17 

Frank, J., King, A., & Raine, D. 1992, Accretion Power in 
Astrophysics. (Cambridge University Press) 

Guilloteau, S. 2005, private communication 

Hubickyj, O. 2001, private communication 

Ida, S. & Lin, D. N. C. 2004a, ApJ, 604, 388 

Ida, S. & Lin, D. N. C. 2004b, ApJ, 616, 567 

Ikoma, M., Nakazawa, K., & Emori, H. 2000, ApJ, 537, 1013 

Kornet, K., Bodenheimer, P., Rozyczka, M., & Stepinski, T. F. 
2005, A&A, 430, 1133 

Kornet, K, Bodenheimer, P., & Rozyczka, M. 2002, A&A, 
396, 977 

Kornet, K, Rozyczka, M., & Stepinski, T. F. 2004, A&A, 417, 
151 

Kornet, K, Stepinski, T. F, & Rozyczka, M. 2001, A&A, 378, 
180 

Laughlin, G, Bodenheimer, P., & Adams, F. C. 2004, ApJ, 612, 
L73 

Laughlin, G. & Rozyczka, M. 1996, ApJ, 456, 279 
Lissauer, J. J. 1993, ARA&A, 31, 129 

Marcy, G. W., Butler, R. P., Fischer, D., et al. 2001, ApJ, 556, 
296 

Marcy, G. W., Butler, R. P., Fischer, D. A., et al. 2002, ApJ, 
581, 1375 

Marcy, G. W., Butler, R. P., Vogt, S. S., Fischer, D., & Lissauer, 
J. J. 1998, ApJ, 505, L147 



